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^^ \ Abstract. The GiBUU model, which implements all reaction channels relevant at medium neutrino 

CN ■ energy, is used to investigate the neutrino and antineutrino scattering on iron. Results for integrated 

cross sections are compared with NOMAD and MINOS data. It is shown, that final state interaction 

can noticeably change the spectra of the outgoing hadrons. Predictions for the Minerva experiment 

are made for pion spectra, averaged over NuMI neutrino and antineutrino fluxes. 
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INTRODUCTION 

b 

p^ \ Neutrino and antineutrino scattering on nuclei for neutrino energies above 30 GeV was 

studied in several experiments starting from the 80s (see [1,2] for a review). Based 
on muon detection, integrated cross sections were measured with combined precision 
i^ \ of 2%; they grow linearly with energy, which agrees with the predictions of the quark 

^ ■ parton model. Double differential cross sections with respect to muon variables were 

O^ ■ also measured, allowing one to extract the nucleon structure functions. 

'^ \ At lower energies the neutrino reactions are not so easy to model, because of the 

C^ \ overlapping contributions from QE scattering, resonance production and background 

processes. Another difficulty is in choosing an appropriate way to describe the onset of 
the DIS processes, the so-called Shallow Inelastic Scattering. This requires complex 
approaches that take all of the relevant channels into account. The lack of data at 
intermediate energies has long been an obstacle for a serious test of such approaches. 
Up-to-date and coming experimental results are changing the situation. Recently NO- 
5^ \ MAD reported the inclusive neutrino cross section for E^ > 4.6 GeV with an accuracy 

of at least 4% [3]. MINOS reported both neutrino (for Ev > 3.48 GeV) and antineutrino 
(for Ey > 6.07 GeV) cross sections with a comparable precision [4]. The Minerva ex- 
periment intends to perform measurements on Plastic (CH), Iron, Lead, Carbon, Water 
and liquid Helium targets in the NuMI beam, which would directly allow to compare 
nuclear effects on various nuclei. Besides muon detection, this experiment will also be 
able to resolve various final states by identifying the tracks of the outgoing hadrons. 

Here we use the GiBUU transport model [5, 6] to study neutrino and antineutrino 
scattering on iron. Our results are compared with the recent MINOS and NOMAD 
data. Predictions are also made for the spectra of the outgoing particles, which can be 
measured by the Minerva experiment. The calculations are done without any fine tuning 
to the data covered here with the default parameters as used in the GiBUU framework. 
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GIBUU TRANSPORT MODEL. VARIOUS NEUTRINO 
PROCESSES AND PROBLEM OF DOUBLE COUNTING 

The GiBUU model was initially developed as a transport model for nucleon-, nucleus-, 
pion-, and electron- induced reactions from some hundreds MeV up to tens of GeV. 
Several years ago neutrino-induced interactions were also implemented for the energies 
up to a few GeV. Recently the code was extended to describe also the DIS processes 
in neutrino reactions. Thus, we can study all kind of elementary collisions on all kind 
of nuclei within a unified framework. The model is based on well-founded theoretical 
ingredients and has been tested against various nuclear reactions. For a detailed review 
of the GiBUU model see [6]. 

GiBUU describes all processes relevant at medium energies. Our approach to 
quasielastic (QE) scattering, resonance (RES) production and background (BG) pro- 
cesses is described in [7, 8]. The DIS scattering is included as Pythia simulation. 

In the region of the shallow inelastic scattering, that is at moderate invariant masses, 
1.6 GeV <W < 2.0 GeV, there is a potential problem of double counting. Here the same 
physical events can be considered as originating from decays of high mass baryonic 
resonances or from DIS. In the GiBUU code we use the ansatz, that both RES, BG and 
DIS processes contribute in this region. The RES and BG contributions are smoothly 
switched off above V7res-2 and DIS contribution is switched on above H^dis-i^ 

,„,x Wres-2-W W-Wdis-i ... 

0-RES = 0-REs(VV)— — 0-DIS = 0-DIs(W)— — (1) 

VVrES-2 — VVRES-l VVDIS-2 — V^DIS-l 

The default parameters are: W^res-i = 2.0 GeV, V7res-2 = 2.05 GeV, H^dis-i = 
1.6 GeV, and W^dis 2 = 1-65 GeV. With this choice, the DIS events become noticeable 
at neutrino energies around 3 GeV. This choice is mainly motivated by the comparison 
with the electroproduction data [9]. 

In essence, the DIS processes below V7res-2 account for the resonances whose elec- 
tromagnetic properties are not known and for the non-resonant processes giving a few 
mesons in the final state beyond the one pion background. 

INTEGRATED CROSS SECTIONS. 

As already mentioned, the DIS cross section grows linearly with energy. So at high neu- 
trino energies the data are conveniently presented as cross section per energy Otot/Ey. 
The world average values are given [10] for isoscalar target: Otot/Ev{v) = 0.667 ± 
0.014-10-38 ^j^2^ Qgy for neutrinos and Otot/Ey{v) = 0.334±0.008 ■ lO^^s ^^2^ Q^y 
for antineutrinos. 

Fig. 1 shows the isoscalar cross section together with the free cross section for the iron 
composition, that is for the 26 protons and 30 neutrons (46.4% protons, 53.6% neutrons). 
Recent MINOS data for iron target as well as NOMAD data for an isoscalar target are 
also shown. 

This figure shows that the decreasing slope of our curves for the neutrino cross section 
is in agreement with that of the data. This slope is not taken into account in deriving 



the world-average value, where it was assumed to be negligible. Notice, that such a 
comparison with the data is meaningful, only if nuclear corrections are neglected. 
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FIGURE 1. (Color online) Total elementary cross section per nucleon for neutrino (left) and antineu- 
trino (right) induced reactions. The world average values [10] and their errorbands are also shown as 
long-dashed green lines. Experimental data are from [3, 4]. 
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FIGURE 2. (Color onUne) Total cross section per nucleon for neutrino (left) and antineutrino (right) 
induced reactions on iron target. Various prescriptions to include nuclear effects in DIS are compared. 
MINOS data (full circles) are from [3, 4], IHEP-JINR data (full diamonds) are from [11]. 

The actual value of such corrections for neutrino reactions is not known so far, because 
of both experimental inaccuracies and difficulties in the theoretical description. On one 
hand, nuclear parton distributions, based on electromagnetic scattering data and intended 
for description of both charged lepton and neutrino reactions, were introduced. For a 
review and a list of recent parametrization see, for example, [12]. On the other hand, 
recent investigation [13, 14] showed, that in neutrino reactions nuclear corrections to 
parton distributions are at the same level as for electrons, but have a very different 
dependence on the Bjorken x variable. The topic remains controversial, with the hope 
that future precise IVIinerva results on various targets will clarify the situation. 

As we already mentioned, the GiBUU code uses PYTHIA for the simulation of DIS 
processes. However, the PYTHIA code was designed for elementary reactions. In the 
GiBUU simulation the neutrino interacts with one initial nucleon, bound in the hadronic 
potential and having nonzero Fermi momentum. In order to be able to use the PYTHIA 
event simulator, we have to provide some quasi-free kinematics as inputs to PYTHIA. 
Various prescriptions to do this (for details see [6]) result in a 5 — 7% difference in the 



results. The corresponding cross sections (denoted as "F-NO", "F-CM", "F-THRE") are 
shown in Fig 2. Nuclear parton distribution functions from [15] are also implemented as 
one of the options to use (to avoid double counting, nuclear potential and Fermi motion 
in such calculations are switched off). The result ("nuclearPDF") as well as the free 
cross section for iron composition ("26p+30n") are also shown in Fig 2. At the moment 
we consider the various prescriptions mentioned above as intrinsic uncertainty of the 
GiBUU code, reflecting the lack of our understanding the nuclear effects. No other event 
generator, as far as we know, accounts for nuclear corrections in high-energy neutrino 
reactions. 

Fig. 2 shows, that our calculations are in a good agreement with the recent neutrino 
data, which are also consistent with each other. For antineutrino the agreement is good 
for Ey > 25 GeV. For lower energies our curve is above the recent MINOS data, but 
below the IHEP-JINR results [11]. The overall agreement of our calculations with the 
data is therefore better than the agreement of the data with each other. 

FINAL STATE INTERACTIONS AND CHANGE OF THE FINAL 

HADRONIC SPECTRA. 
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FIGURE 3. (Color online) Pion kinetic energy distributions for neutrino and antineutrino induced 
reactions for 1-pion production (one pion of a given charge and no other pions are produced) "MULTI" 
at least 1 pion of a given charge and any number of pion of other charges are produced". Calculations are 
for the NuMI low-energy mode neutrino/antineutrino fluxes. 

After being produced in the initial interaction act, outgoing hadrons propagate 



throughout the nucleus. In GiBUU this process of final state interactions (FSI) is mod- 
eled by solving the semi-classical Boltzmann-Uehling-Uhlenbeck equation. It describes 
the dynamical evolution of the phase space density for each particle species under the 
influence of the mean field potential, introduced in the description of the initial nucleus 
state. Equations for various particle species are coupled through this mean field and also 
through the collision term. This term explicitly accounts for changes in the phase space 
density caused by elastic and inelastic collisions between particles. 

FSI decrease the cross sections as well as significantly modify the shapes of the final 
particle spectra. Such change was seen, for example, in photo-pion production [16] and 
is described by the GiBUU with a good accuracy. A similar change should be observed 
in neutrino reactions. 

Fig. 3 shows the n^ (upper panels), n^ (middle panels) and K~ (lower panels) spectra 
for neutrino (left panels) and antineutrino (right panels) NuMI fluxes for 1-pion events. 
The green dashed lines show the kinetic energy (Tj,;) distributions without FSI, i.e. of 
pions produced in the initial neutrino vertex. The solid red lines show the distributions 
after FSI, i.e. of pions that made it out of the nucleus. Such spectra can also be calculated 
for any other predefined final state and should be measurable in Minerva experiment. 

For dominant channels (;r+ production for neutrino reactions and ;r in antineutrino 
ones), the FSI decrease the cross section at 7,^ > 0.2 GeV. This is mainly explained 
by pion absorption through nN — > A following by AN -^ NN. Pion elastic scattering 
in the FSI also decreases the pion energy, thus depleting the spectra at higher energies 
and accumulating strength at lower energies. Thus, an increase of the cross sections is 
observed at T;,: < 0. 15 GeV; altogether this leads to a significant change of the shape of 
the spectra. 

Scattering can also lead to pion charge exchange. For neutrino-induced reactions, the 
n^n -^ n^p scattering in the FSI is the main source of side-feeding for the n^ channel, 
leading to a noticeable increase of the n^ cross section at low Tj,;. The inverse feeding 
is suppressed, because less n^ than ;r+ are produced at the initial vertex. The same 
mechanism of side feeding from dominant to sub-dominant channel through n^ p -^ K^n 
is working for antineutrino induced reactions. 

For the least dominant channel {tz^ production for neutrino reactions and ;r+ in 
antineutrino ones), the FSI (in particularly side feeding) represent the main source of 
the events observed; thus a dramatic FSI effect. 

Fig. 4 shows the origin of the pions (that is the initial vertex, at which the pion 
was produced) in the dominant channels for various final states. 1-pion production, for 
example, receives a major contribution comes from Delta resonance production and its 
following decay. For other final states with more pion DIS dominates, but the Delta is 
still visible. The contribution from the QE vertex is very small but nonzero, the outgoing 
pion in this case can only be produced during the FSI, for example, due to the NN — > NA 
scattering followed by A ^ Nn. 

Similar plots can be obtained from the GiBUU simulation for any other outgoing 
particles (protons, kaons, eta) for any predefined final state and compared with the 
coming Minerva results. 
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FIGURE 4. (Color online) Pion kinetic energy distributions for 7r+ production in neutrino reaction and 
n^ production in antineutrino reaction (both are dominant channels), showing various contribution to a 
given final state. "1-pion": one pion of a given charge and no other pions are produced in the final state 
after FSI; "IX": one pion of a given charge and any number of other pions; "2X": two pions of a given 
charge and any number of other pions; "MULTI": at least 1 pion of a given charge and any number of 
other pions. Calculations are for NuMI low-energy mode neutrino/antineutrino fluxes. 
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